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ABSTRACT 

The microwave "haze" was first discovered with the initial release of the full sky data from the 
Wilkinson Microwave Anisotropy Probe. It is diffuse emission towards the center of our Galaxy 
with spectral behavior that makes it difficult to categorize as any of the previously known emission 
mechanisms at those wavelengths. With now seven years of WMAP data publicly available, we have 
learned much about the nature of the haze, and with the release of data from the Fermi Gamma- 
Ray Space Telescope and the discovery of the gamma-ray haze/bubbles, we have had a spectacular 
confirmation of its existence at other wavelengths. As the WMAP mission winds down and the Planck 
mission prepares to release data, I take a last look at what WMAP has to tell us about the origin of 
this unique Galactic feature. Much like the gamma-rays, the microwave haze/bubbles is elongated in 
latitude with respect to longitude by a factor of roughly two, and at high latitudes, the microwave 
emission cuts off sharply above ~35 degrees (compared to ~50 degrees in the gammas). The hard 
spectrum of electrons required to generate the microwave synchrotron is consistent with that required 
to generate the gamma-ray emission via inverse Compton scattering, though it is likely that these 
signals result from distinct regions of the spectrum (~10 GeV for the microwaves, ~1 TeV for the 
gammas). While there is no evidence for significant haze polarization in the 7- year WMAP data, I 
demonstrate explicitly that it is unlikely such a signal would be detectable above the noise. 

Subject headings: Galaxy: center — ISM: structure — ISM: bubbles — Radio continuum: ISM 



1. INTRODUCTION 

With the initial release of the full sky data by 
the Wilkinson Micr owave Anisotropy Probe (WMAP, 
iBennett et al.ll2QQ3bl ) came the discovery that there was 
anomalously hard spectrum microwave emission towards 
the center of our Galaxy (Fin kbeinerl l2QQ4at ). The 
method for uncovering this WMAP "haze" as it was 
termed at the time was very simple: there were known 
or anticipated microwave emission mechanisms from the 
Galaxy and it was believed that templates (maps at 
other wavelengths) morphologically tracing these fore- 
grounds could be used to clean Galactic emission from 
the maps for CMB analysis and to study the foregrounds 
in their own right ([Bennett et al.ll2003af ) using straight- 
forward linear regression techniques. These templates 
were known to be insufficient down to the measurement 
noise in the WMAP maps, yet the template fitting proce- 
dures worked remar kably well for both the cosmological 
([Spergel et al.ll2QQ3h and Galactic science purposes. 

Ho wever, while removing; >95% of the variance in the 
data (|Finkbeinerll2QQ4aHDobler fc Finkbeinerl l2QQ8a) the 
templates failed most spectacularly within ~ 30° of the 
Galactic center (GC). After removing the emission cor- 
related with the templates, the residual resembled a 
"hazy" blob with a spectrum that was too soft to be free- 
free emission and too hard to be synchrotron emission 
from electrons generated by supernova (SN) remnants 
(jDobler fc Finkbeinerl l2008al hereafter DF08). This mi- 
crowave haze is somewhat elongated in Galactic latitude 
with an intensity falling off with distance from the GC. 
A careful measurement of the spectrum and radial de- 
pendence of the haze by DF08 led to ample speculation 
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about the origin of the h aze electrons: from st ar forma- 
tion and SN acceleration ([Biermann et al. 2010) to signa- 
tures of particle dark matter annihilation in o ur Galactic 
halo (iHooper et aLll2QQ7t ICholis et afll2Q09allbh to claims 
that the feature did not exist at all ([Mertsch fc Sarkarl 

[2oToh . 

The last of these claims was laid to rest with the 
first year data release of the Fermi Gamma-Ray Space 
Telescope w hich showed a simila r structure at gamma- 
ray energies Dob ler et al.l (|2010l ). The implication was 
that the same electrons that generate microwaves at 
WMAP wavelengths generate gammas via inverse Comp- 
ton (IC) scatt ering; that are obser ved by Fermi. It was 
also found by iDobler et al.l ([20 10) that, while the ex- 
tent of the microwaves is <~ ±30°, the gamma-ray haze 
extended to a full ±50° above and below the Galac- 
tic plane, making it the second largest structure in our 
Galaxy after the Galactic disk (assuming that it is cen- 
tered roughly on the GC) with a top-to-bottom length 
of 20 kp c. Furthermore , analysis of 1.6 years of Fermi 
data by ISu et al.l (|201Q[ ) suggests that the gamma-ray 
emission also has sharp edges and it was renamed the 
"bubbles". Given the established correspondence be- 
tween the WMAP "haze" and the Fermi "haze/bubbles", 
I will refer to the microwave emission as the WMAP 
"haze/bubbles" or "haze" interchangeably throughout 
this paper. 3 

We now have the benefit of six additional years of mi- 
crowave data from WMAP and clues from the Fermi 
gammas. Given that the release of the Planck data is 

3 Historically, there has been some ambiguity to the nomencla- 
ture and to what the term "haze" refers. Given the observations, 
the clearest definition would be that the haze is a population of 
anomalously hard spectrum cosmic-rays towards the GC and the 
microwave and gamma-ray haze/bubbles are observed signals gen- 
erated by those cosmic-rays. 
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scheduled for late 2012, I will present a description of 
what we know about the haze/bubbles at the end of the 
WMAP era and the beginning of the Planck era. In par- 
ticular, in $2] I will describe the methods used to uncover 
the haze in the first place and comment on other com- 
ponent separation techniques. In §3] I present the 7- year 
WMAP haze/bubbles spectrum and morphology and de- 
scribe new findings at high latitude. In ^4] I address the 
state of the haze in the WMAP polarization data and in 
$5] I discuss the merits and demerits of the various pro- 
posals for the origin of the haze/bubbles electrons before 
summarizing in §6\ 

2. COMPONENT SEPARATION 

Component separation in the context of CMB fore- 
grounds is a broad term for identifying and remov- 
ing (separating) Galactic emission mechanisms from the 
data. Template fitting is perhaps the most basic of 
the available techniques but has proven to be remark- 
ably powerful despite its simplicity. The fundamental 
equation to be solved is just a linear fit of templates to 
the data: w = Pa where w is a map of the WMAP 
data, P is a matrix of template maps, and a is a vec- 
tor of amplitudes. The least-squares solution to this 
equation (after taking the noise map n into account) 
is a = (P T n _1 P) _1 (P T n _1 w). This formalism easily 
accommodates partial sky coverage and arbitrary num- 
bers of templates. Of course, for the purpose of study- 
ing Galactic foregrounds, one of the templates must be 
a CMB estimate with an i nput CMB spectrum, and as 
shown in DF08 (a s well as iDobler fc Finkbeinerl l2QQ8ri 
iDobler et al.ll2009h . this introduces a bias in the derived 
foreground spectra because no CMB estimate is ever 
completely clean of the very foregrounds to be measured. 

When applying template fitting to the WMAP 7- year 
temperature data using the templates and partial sky re- 
gions defined in DF08 (to take into account foreground 
spectral variation with position), the 7- year WMAP 
haze /bubbles emission is shown in Figure [TJ However, 
I have made two adjustments to the templates used in 
the fit compared to DF08. First, for the haze template, 
I use a bivariate Gaussian of scale length = 15° and 
a b = 25°. Second, the haze residual in DF08 included 
a significant "disky" component near the mask likely as- 
sociated with energy dependent propagation lengths of 
ele ctrons accelerated in the Galactic disk as discussed 
in iMertsch fc Sarkarl ([201(f). This is a failing of using 
the Haslam 408 MHz map ([Haslam et al.l ll982) as a soft 
synchrotron template. Additionally, there is the possi- 
bility that a population of pulsars in the Galactic disk 
could potentially produce harder spectrum cosmic-ray 
electrons. 4 Whatever its origin, this disky component 
likely does not represent true haze emission and as such 
I include a Gaussian disk template with scale lengths 
<j£ — 20° and = 5° . The mask used to fit the data ex- 
cludes regions where the dust extinction at Ka is greater 
than 1 magnitude as well as all point sources in both the 
WMAP and Planck ERCSC (30 GHz to 143 GHz) cata- 
logs, the LMC, SMC, M31, Orion-Barnard's Loop, NGC 

4 However, the efficiency factor for generating electrons with 
pulsars is prese ntly not constrained (and could be as low as zero), 
while Mertsch & Sarkar (2010) showed that the energy dependent 
propagation effects are of roughly the right order of magnitude 
given the results in 



5128, and C-Oph. 

As seen in Figure [TJ outside of the haze region, the 
residuals are remarkably small, indicating that the tem- 
plates used in the analysis do represent reasonable trac- 
ers of the actual emission. 5 Most importantly, several 
authors have found that emission from rapidly rotating 
dust grains (which have small dipole moments and thus 
produce spinning dipole radiation) represen ts >30% of 
the d iffuse emission at K-band (see DF08 and lFinkbeinerl 
2004a). This spinning dust is perhaps the most uncertain 
foreground since it is not well known how its emissiv- 
ity scales with total dust column density (inferred from 
100/im for example). 

In addition to the simple template fitting described 
above, more involved component separation techniques 
have been developed and applied to the WMAP data. 
The two most commonly used are the Maximum En- 
tropy Method ([Bennett et al.l l2003af) and Baye sian in- 
ference via Gib os sampling (|Eriksen et al.ll2006h . While 
these methods are excellent for cleaning the CMB of fore- 
grounds and estimating the amplitude of emission with 
known morphology (such as the cosmic dipole), they are 
not optimal for identifying new components. For ex- 
ample, although it has been reported that these tech- 
niques do not " s ee" th e WMAP haze/bubbles, in fact 
iPietrobon et al.l (|201ll ) show that the haze emission is 
typically swept into a "low frequency" component that 
is a conglomerate of spinning dust, soft synchrotron, hard 
synchrotron, and free-free all forced to obey a single 
power law spectrum. 

3. SPECTRUM AND MORPHOLOGY 

3.1. The Haze at Low Latitudes 

The residuals in Figure [T] and the spectrum plots in 
Figure [2] show the two important features of the mi- 
crowave haze/bubbles. First, the spectrum is very hard, 
Th oc v$ h with Ph ~ —2.5 which is harder than can 
be obtained from particle acceleration in SN shocks af- 
ter taking into account cosmic-ray diffusion and energy 
losses. However, it is significantly softer than the /3p = 
—2.15 required if the emission were free- free. Second, the 
haze/bubbles are extended and elongated in latitude. In 
detail however, while the haze is roughly centered on the 
GC, its morphology is somewhat different in the north 
and south. In the south, and above latitude b = —30°, 
the morphology of the WMAP haze/bubbles is strikingly 
similar to the Fermi haze/bubbles (right panel of Figure 
[T|). In the north, the situation is much more complicated 
by the significant dusty regions of the Ophiuchus com- 
plex. In fact, the direction of the north Galactic center 
(below b = 30°) likely contains either significant spin- 
ning dust emission that is not wel l appr oximated by our 
dust template (the Schleg eTet al.l (Il998[) dus t map eval- 
uated at 94 GHz bv iFinkbeiner et al.l (|l999h ) or highly 
variable gas temperatures making our free-free template 

5 The one exception is the Gum Nebula at t ~ —1 00°; there are 
two re asons for this. First, as pointed out in DF08 and Doble r~et al.l 
(2009), the gas temperatures and spinning dust spectra are vary- 
ing very rapidly with position in this region. Second, the Gum 
Nebula is one of the brightest Galactic regions in the raw data and 
so these residuals actually represent a small fraction of the total 
emission. Also, it is important to note that the Gum Nebula is rel- 
atively nearby and as such occupies a much smaller volume than 
the haze/bubbles. 




Ka-band 7-year WMAP haze/bubbles Q-band 7-year WMAP haze/bubbles 




Fig. 1. — The haze/bubbles in both microwaves and gamma-rays (23 GHz, 33 GHz, 41 GHz, and 2-5 GeV, counter clockwise from top 
left). In the microwaves, templates have been used to regress out emission from the CMB, thermal an d spinning dust, fre e-free, and soft 
synchrotron. In the gammas, the official Fermi diffuse model has been subtracted from the data (see Dobler et al. 2010). In all bands, 
the haze is seen to be elongated in latitude by a factor of roughly two reaching ±50° in the gammas and ~ ±35° in the microwaves. The 
microwaves are stretched with a z/ 2 ' 5 scaling, which yields roughly equal brightnes s from K to Q ban d, indicat i ng an electron spectrum 
of dN/dE e oc E~ 2 , broadly consistent with the gamma-ray spectrum from IDobler et al. ( 2010) and Su et al. ( 2010|). The gamma-ray 
haze/bubbles seem to have a sharp edge near |6| ~ 50° while the microwaves seem to fall off quickly for |6| > 30° (particularly in the south). 




Fig. 2. — Left and middle: scatter plots (drawn with contours) for the microwave haze/bubbles residuals in Figure^ (dashed lines) and 
for the total synchrotron (haze plus soft synchrotron; solid lines) in the region \l\ < 25°, —35° < b < —10°. The best fit power law for 
this region shows that the haze emission is significantly harder than the soft synchrotron from 23 GHz to 41 GHz. The total emission in 
the region for both cases is shown in the right panel. The haze emission is clearly harder given the noise in the data, is consistent with a 
power law of roughly z/ -2-5 , and represents approximately 33% of the total synchrotron emission at K-band. 



(the iFinkbeinerl (j2QQ3h Ha map) unreliable. In reality, 
both effects are likely operating simultaneously. 

This additional complexity in the north is the reason 
that I concentrate on the south when estimating the 
spectrum of the haze in Figure El The derived spec- 
trum of Ph ~ —2.5 implies an electron spectrum of 
dN[dE e oc E~ 2 which ha s been shown by IDobler et al.l 
(|2010D and lSu et aD (|2010f ) to be broadly consistent with 
the spectrum of the Fermi gamma-ray haze/bubbles. 
This, taken together with the morphological correlation 
at low latitudes, provides the strongest evidence that 
this is in fact the same phenomenon observed at mul- 
tiple wavelengths. While it is true that DF08 showed 



that the bias induced by presubtracting a CMB esti- 
mate can be quite large, for a component with a spec- 
trum T oc u~ 2 - b the CMB5 weights for the WMAP 7 
data (C = [0.246,-0.736,-0.0685,0.263,1.295] in ther- 
modynamic AT for the five WMAP bands) would lead 
to a measured spectrum ex i>~ 2 - 53 when comparing K- 
to Ka-band and oc ^~ 2 - 54 when comparing K- to Q- 
band. In other words, the bias for a foreground with a 
haze/bubbles-like spectrum would be rather small given 
the DF08 CMB5 estimate and so (3h = -2.5 is likely very 
close to the true spectrum. Furthermore, a component 
with a free-free spectrum would be inferred to have a 
spectrum of oc ^ _2154 providing more evidence that the 
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Fig. 3. — Top and middle: intensity in CMB IC scattered gam- 
mas and microwave synchrotron as a function of the electron energy 
for a population of electrons with spectrum dN/dE e oc E~ 2 illus- 
trating that the WMAP haze/bubbles and Fermi haze/bubbles are 
generated by different ranges of the electron spectrum. In particu- 
lar, if the gammas at ±10 kpc are generated by scattering of CMB 
photons, this requires ~TeV electrons at those distances, while the 
microwaves are generated by electrons with energies ~10 GeV. The 
higher energy electrons used to generate the gammas would create 
a synchrotron signal that peaks at very high frequencies (bottom), 
though it would be overshadowed by thermal dust emission. 



microwave haze/bubbles is not due to free-free emission. 

However, it is important to bear in mind that the ma- 
jority of the haze/bubble gammas observed by Fermi at 
high latitudes likely come from E e > 100 GeV electrons 
scattering CMB photons, while the microwaves are gen- 
erated by electrons with energies 1 GeV < E e < 100 
GeV. This is illustrated in Figure [3] which shows the to- 
tal emissivity 



•tot 
Jv 



J E e j v (E e ,B)^d\nE e 



(1) 



where j v is the emissivity from a single electron at energy 
E e in a magnetic field B, as a function of electron en- 
ergy for a population of electrons with the haze/bubbles 
spectrum in a 5/iG magnetic field. It is also worthwhile 
to note that, as shown in the figure, if the electron en- 
ergy extends up to order TeV energies as required by 
the gammas, then the synchrotron emissivity per elec- 
tron for E e > 100 GeV actually increases with increas- 
ing frequency. However, the falling electron spectrum 
and emission from thermal dust (which rises much faster 
with frequency) render this high frequency synchrotron 
undetectable. 

3.2. The Haze at High Latitudes 

An aspect of the haze/bubbles that has not been fully 
discussed is the behavior of the microwave emission at 
high latitudes, where the gamma-ray signal is most un- 
ambiguous and where the microwaves appear to fall off in 
intensity. Figure |4] shows the microwave and gamma-ray 
haze side-by-side, with a large ±28° mask and smoothed 
to a common 2° beam. It is immediately clear from the 
figure that there most certainly is haze emission above 
\b\ = 28°, but its morphology is somewhat different than 
the gammas. 

Concentrating again mostly on the south, the gam- 
mas fill a bubble with a sharp edge while the microwaves 
seem to fall off dramatically for b < —35° as shown in the 
lower left panel of Figure |U The intensity as a function 
of latitude within the bubble (defined by the gamma-ray 
emission), does not quite have an "edge" in the same 
sense as the gammas, but the fall off is quite rapid in- 
dicating either a rapid decrease in the number of elec- 
trons or the strength of the magnetic field. Given the 
gammas at latitudes —50° < b < —35°, the latter seems 
more likely. Furthermore, given that neither synchrotron 
maps at lower frequencies (e.g., the 408 MHz map) nor 
WMAP polarization show a sharp drop in synchrotron 
emission below b = —35°, this seems to indicate that the 
dominant field component producing the synchrotron is 
the field within the haze/bubbles itself rather than the 
Galactic magnetic field. Assuming that the haze/bubbles 
is located at the GC, then the scale height of this mag- 
netic field drop off is ~6 kpc. 

The "sharpness" of the microwave edge is depicted in 
Figure [5] which shows the total intensity as a function of 
both latitude and longitude in the south. Interestingly, 
for latitudes —35° < b < —6°, the brightness is roughly 
flat with distance. There is some slow decrease, but the 
drop below b = —35° is striking, and this flatness is rem- 
iniscent of the gamma-ray emission — though again, the 
edge is not as sharp and is at lower absolute latitude. 
This latitudinal profile is qualitatively different from that 
presented in DF08 in which the haze brightness increases 
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Fig. 4. — The Fermi (top left) and WMAP (top right) haze/bubbles at high latitudes masking emission for |6| < 28°. The correlation 
between the two is clear, though the microwave "edge" is at a lower latitude than the gamma-ray edge, indicating that the "bubble" is not 
completely filled by microwaves. The dashed white line at £ = —5° illustrates that, while the haze/bubbles is centered on the GC at low 
latitudes, it is offset by a few degrees above |6| ~ 28°. The sharp drop in microwaves in the Galactic south (bottom left) indicates either 
a sharp decrease in the total number of electrons or a sharp drop in the magnetic field. The former seems unlikely given the presence of 
gammas down to b ~ —50°. Bottom right the same stretch, but with contours showing T an tenna = 10 -2 and 10 -3 mK indicating the very 
sharp transition from haze/bubbles to noise below b ~ —30° (see also Figure 0). 



sharply towards the GC above b = —10°. The difference 
is that the inclusion of the extra disk template (to ac- 
count for the excess disky emission near the mask) has 
removed most of this power making the profile "flatter" 
with latitude. As a function of longitude, again the haze 
brightness is roughly constant within about \£\ < 15° 
outside of which the fall off is quite rapid. The figure 
also shows that for latitudes —50° < £ < —40°, there is 
no significant microwave haze/bubbles emission, in con- 
trast to the gammas. 

In the north, again the situation is somewhat compli- 
cated since there is clearly emission co-located in both 
the microwaves and gammas, but there is also significant 
contamination by spinning dust in the microwaves and 
7T° decay photons in the gammas. Nevertheless, it is clear 
from the morphology (see Figure H]) that, while the lower 
latitude microwaves are centered on the GC, the center 
of the structure at \b\ = 28° is roughly £ = —5°, an offset 



which causes some difficulty for formation scenarios. 

4. POLARIZATION 

Given the spectrum derived in Figure [2j the haze emis- 
sion certainly seems to be consistent with hard spec- 
trum synchrotron. Also given that synchrotron with 
/3 ~ — 2.5 can be up to 70% polarized in a completely or- 
dered magnetic field, this raises the question whether the 
haze appears in WMAP polarization data. Several au- 
thors have addressed this issue, most notably [Gold et ahl 
(2011) who claim no evidence for a significant hard com- 
ponent from WMAP 7- year data. However, this claim 
suffers from two major drawb acks. First, as shown in 
iMiville-Deschenes et all (|2QQ8f ) and argued in DF08, even 
a small turbulent component in the magnetic field can 
reduce the polarization amplitude when projecting along 
the line of sight. In fact, this line of sight depolariza- 
tion is clearly observed in the WMAP data as illustrated 
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Fig. 5. — The microwave haze/bubbles profile as a function of latitude (left) and longitude (right). Inside of the haze region, the brightness 
is roughly flat (or at most, decreases slowly with latitude) and then drops sharply for latitudes below b « —35° and for longitudes \£\ > 15°. 
This "edge" in the microwaves is not as sharp as in the gammas, and the dotted line in the right panel indicates that there is little to no 
microwave emission for latitudes —50° < b < —40°, in contrast with the gamma-ray measurements, indicating a relatively quick drop in 
the magnetic field at b ~ —35°. The dotted line in the left panel shows the results without including the bivariate Gaussian disk template. 
In that case, the profile rises dramatically towards the GC, but this emission is not likely to be associated with the haze (see 



in Figure [6] which shows that (outside of the mask used 
in the fits), almost all of the emission is due to Loop I 
which is a local feature. The disk synchrotron is sim- 
ply not present in the polarization data though it is seen 
in the total intensity data, indicating that most of the 
disk emission is line-of-sight depolarized by the varying 
orientation of the Galactic magnetic field. 

As Figure [6] illustrates, the second shortcoming is that, 
even if the polarization fraction of the haze were constant 
with respect to the polarization fraction of prominent 
features like the Loop I SN remnant, the noise in the 
polarization data is prohibitively large to detect a harder 
synchrotron component. The argument goes as follows, 
forming the polarization residual 

, x-3.3 

Rp=Pk,-PkX r-^\ , (2) 

where Pi and v\ are the WMAP polarization map and 
frequency at band i respectively, produces a map consis- 
tent with noise as shown in Figure [6] That is, all of the 
emission (synchrotron in the case of polarization) out- 
side of the mask is consistent with having a single power 
law spectrum j3s = —3.3. However, if we make the same 
residual using just the synchrotron intensity Si, 

R I = S Ka -S K x(^j , (3) 

where Si is the template extracted synchrotron (i.e., the 
raw data minus the CMB, thermal plus spinning dust, 
and free- free components), Figure [6] shows that this resid- 
ual also shows no evidence of the haze. However, we 
know that the haze is there and that its spectrum is 
harder than elsewhere in the Galaxy from Figures [TJH] 
and Figure 7 of DF08. The conclusion is that, given the 



short lever arm of 23-33 GHz, even the temperature data 
are too noisy to extract a slightly harder spectrum com- 
ponent in this manner, let alone the polarization data. 

5. ORIGIN SCENARIOS 

Since the haze/bubbles was f irst discovered in the mi- 
crowaves by iFinkbeinerl ([2004aj ) , an d especially since th e 
recent discovery in the gammas bv lDobler et all (|2010h . 
there has been significant effort devoted to theorizing 
about the origin of this hard spectrum population of 
electrons. Each scenario has associated pros and cons, 
successfully producing some aspects of the emission and 
failing to produce others. The following is a list of the 
leading possibilities put forth in the literature. 

Galactic wind: The Galactic win d scenario 

(jCrocker et al.l 1201 It iCrocker fc Aharonianl l2011ft sug- 
gests that cosmic-rays are accelerated to ~TeV energies, 
that the gammas are the result of collisions of cosmic-ray 
protons with a very underdense ISM in the bubbles pro- 
ducing 7T° decay emission, and that the microwaves rep- 
resent the synchrotron from secondary electrons. This 
model can reproduce the hard spectrum WMAP signal 
and, because the proton spectrum is also very hard, the 
hard spectrum gamma-rays as well. The primary fail- 
ing of this model is that it requires extended injection 
of cosmic-rays over several billions of years, making the 
sharp edge of the gammas and microwaves very difficult 
to maintain over those long time scales. Additionally, 
even in the event where the ISM is "saturated" with pro- 
tons as they describe, the bubbles will produce constant 
volume emissivity leading to limb darkening of the gam- 
mas contrary to observ ations. Finally, as in observed 
winds in other galaxies (Veill eux et al.ll2005l ), we would 
expect to see a significant cool component producing Ha 
emission, but no Ha associated with the haze/bubbles 





Fig. 6. — The WMAP polarization (left) and total synchrotron intensity (right) at K- and Ka-band (top two rows). Although some 
features are similar (namely Loop I in the North), there are significant regions of synchrotron intensity that do not appear in polarization 
(e.g., the Galactic plane) indicating line-of-sight projection effects through a turbulent magnetic field that reduce the total polarized signal. 
Scaling Pk by a single (3 = —3.3 power law across the whole sky and then subtracting from Pxa yields residuals consistent with noise 
(bottom left panel) potentially indicating that the haze does not appear in polarization. However, a similar exercise with the total intensity 
also yields similar results (bottom right panel) despite the fact that the haze/bubbles are there and have a different spectrum (see Figures ^ 
and[2j). The implication is that the 23-33 GHz lever arm is too small and the noise in the polarization data too large to positively identify 
or rule out a polarized haze/bubbles component. 



has yet been observed (DF08: iSu et al.ll2010D . 

Starburst/SNe: In this scenario a significant outflow 
is ge nerated by star formation and/or SNe in the GC 
(e.g., iBiermann et al.lfeOlOh . However, this model would 
not only likely lead to associated Ha, but it is not clear 
how the uniform intensity gammas are produced or the 
sharp edges. Furthermore, the timescale for the cosmic- 
rays to reach 5-10 kpc as required by the microwave and 
gamma-ray data implies that the spectrum would be 
rendered too soft to reproduce the haze/bubbles given 
that the energy loss time for ~ 100 GeV electrons is 
r ~ 10 5 — 10 6 yr. Reacceleration of the cosmic-rays is 
a possibility, but what the reacceleration mechanism is 
above a few kpc is unclear. 

Second order Fermi acceleration: Rather than an 
actual "formation" mechanism, this model relies on the 
haze/bubbles to be generated by, for example, a GC jet 
(see below) which is then filled with magnetosonic tur- 
bulen t waves that accelerate particles ([Mertsch fc Sarkarl 
l201lh . At the source of injection, the electron spectrum 
is dN/dE e oc E~ 2 . This scenario suffers from having 
relatively little predictive power since the absolute num- 



ber of accelerated electrons relies on the microphysics of 
the acceleration mechanism (which is not well known) 
and the injection morphology can be set somewhat arbi- 
trarily (though it must be concentrat ed near the outer 
shock ). Furthermore, as found by Mer tsch fc Sarkarl 
(|201ll ). for the required parameters to explain the gam- 
mas, the WMAP haze is underpredicted by an order of 
magnitude. This tension can be alleviated somewhat by 
increasing the magnetic field strength, but this would 
predict microwaves below b = —35° (and indeed a bub- 
ble edge in the microwaves at \b\ ~ 50°) which is not seen 
in the data. 

Active Galactic nucleus: One of the most promis- 
ing scenarios is that there was some accretion event onto 
the central Galactic black hole roughly 1 million year s 
ago that resulted in a GC jet (|Guo fc Mathews! l20Tlh . 
This model is attractive mainly because it is episodic 
(meaning that it can more easily explain the sharp edges 
in the gammas) and it can produce roughly the required 
integrated power in gammas and microwaves. There are 
four primar y concern s with the model however. First, 
as shown by ISu et all (|2010l ) the intensity profile of the 
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gammas is roughly uniform and, at best, the AGN model 
would predict a uniform volume emissivity, again yielding 
limb darkening, especially near the base of the jet gen- 
erated bubble which is very thin. Second, generically, 
shearing instabilities are generated at the edge which 
are not seen in the gamma ray data which has smooth 
edges. However, viscosity withi n the plasma ca n effec- 
tively suppress the instabilities (|Guo et al.ll201l[ ). Addi- 
tionally, magnetic draping may suppress the instabilities 
(and have the added benefit of confining the cosmic-rays 
as required), though it has not yet been shown that the 
required field strengths would not produce synchrotron 
at latitudes \b\ > 35° which is not observed. Relatedly, 
the third concern is that, from jets in other galaxies, we 
see radio lobes - large bubbles of r adio emission at high 
latitudes (e.g.. ISchreier et aDll979h . As shown in Figure 
[U however, the WMAP haze is confined to lower latitudes 
and large radio lobes coincident with the gamma-rays are 
not present. Finally, while the integrated power is suffi- 
cient to explain the haze, the model has the significant 
disadvantage that it has yet to predict a spectrum for the 
accelerated electrons, which is a key characteristic of the 
haze. The softening argument from above still applies 
here and the requirements are that the spectrum 10 kpc 
away from the event have dN/dE e oc E~ 2 . 

Dark matter annihilation: In this model, the dark 
halo of the Milky Way is composed of particles which 
have a self-annihilation cross-section and number den- 
sity sufficient to explain the required injected power. 
This scena r io was first explored in the micr owaves by 
iFinkbeinerl (|20Q4bh and iHooper et all (j20Q7h and then 
expanded to include local cosmic-r ay measurem e nts by 
ICholis et all (|2QQ9bh . Recently, iDobler et all ([201 if ) 
showed that the gamma-ray spectrum, amplitude, and 
morphology can also be reproduced with dark matter an- 
nihilation if the dark halo is prolate and diffusion occurs 
preferentially along ordered field lines in the GC. Here 
the one significant failing of this model is that it does 
not produce sharp edges as seen in the Fermi data (the 
microwave morphology can be completely dominated by 
the magnetic field geometry). Thus, if these sharp edges 
persist with future data, the dark matter annihilation 
only model will be disfavored. 

Of course, it is important to bear in mind that there 
may be multiple mechanisms operating at once (e.g., a 
bubble blown by a jet and filled with electrons by dark 
matter annihilation or a GC wind). In fact, this may 
be the only way to reconcile the very unusual features 
seen in the WMAP and Fermi data, especially given the 
fiat brightness profile of the haze/bubbles which is not 



reproduced in any of the above models. 

6. SUMMARY 

I have presented a last look at the WMAP 
haze/bubbles from the WMAP 7- year data, prior to the 
Planck data release. The haze morphology and spectrum 
are similar to previous analyses with several important 
differences. The inclusion of a disk template in the re- 
gression analysis suggests that the haze profile is some- 
what flatter with latitude than previously thought. Fur- 
thermore, the sensitivity with seven years of data does 
permit a study of the emission at high (southern) lat- 
itudes. Here I find that there is haze emission at lat- 
itudes \b\ > 28° and that there is a relatively sharp 
fall off in intensity for b < —35°, most likely indicat- 
ing a rapid fall off in the magnetic field strength at ~6 
kpc within the haze/bubbles itself. As previously re- 
ported, the spectrum of the WMAP haze is too soft to 
be free-free emission and too hard to be generated by 
SN shocks (the "normal" particle acceleration mecha- 
nism) given diffusion effects. The spectrum is consis- 
tent with that required to produce the Fermi gamma- 
ray haze/bubbles, though I have illustrated that the two 
emissions at high latitudes are likely coming from mostly 
distinct parts of the cosmic-ray spectrum. While the haze 
is not strongly seen in the 7-year WMAP polarization 
data, I have shown that, even if the haze were not de- 
polarized by turbulence in the magnetic field, the noise 
in the WMAP polarization data is sufficient to render 
a hard spectrum polarized component undetectable by 
comparison of K-band and Ka-band. 

Lastly, everything that we have discovered to date 
about the haze/bubbles (hard spectrum synchrotron and 
gammas, sharp edges in the gammas, ±10 kpc extent 
of the electrons, microwaves confined to lower latitudes, 
uniform intensity in gammas, lack of strong associated 
Ha and strong polarization, low energy cutoff in elec- 
tron spectrum, etc.) has made pinpointing a single un- 
derlying origin for the electron population very difficult. 
All of the proposed mechanisms have associated prob- 
lems and if all of the haze/bubbles characteristics persist 
with future data, hybrid formation scenarios will likely 
be required. 
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